Functional Bi2O3/Gd2O3 Silica-Coated Structures for Improvement of Early Age and Radiation Shielding Performance of Cement Pastes

This study presents a new approach towards the production of sol-gel silica-coated Bi2O3/Gd2O3 cement additives towards the improvement of early mechanical performance and radiation attenuation. Two types of silica coatings, which varied in synthesis method and morphology, were used to coat Bi2O3/Gd2O3 structures and evaluated as a cement filler in Portland cement pastes. Isothermal calorimetry studies and early strength evaluations confirmed that both proposed coating types can overcome retarded cement hydration process, attributed to Bi2O3 presence, resulting in improved one day compressive strength by 300% and 251% (depending on coating method) when compared to paste containing pristine Bi2O3 and Gd2O3 particles. Moreover, depending on the type of chosen coating type, various rheological performances of cement pastes can be achieved. Thanks to the proposed combination of materials, both gamma-rays and slow neutron attenuation in cement pastes can be simultaneously improved. The introduction of silica coating resulted in an increment of the gamma-ray and neutron shielding thanks to the increased probability of radiation interaction. Along with the positive early age effects of the synthesized structures, the 28 day mechanical performance of cement pastes was not suppressed, and was found to be comparable to that of the control specimen. As an outcome, silica-coated structures can be successfully used in radiation-shielding cement-based composites, e.g. with demanding early age performances.


Introduction
Concrete is the most popular material used for radiation shielding due to low costs, global availability of its constituents, sufficient durability and its ability to form any shape required [1,2].Moreover, the properties of cement-based composites can be easily altered by the introduction of various heavyweight fine and coarse materials [3,4].Thanks to versatility in designing the mixture composition, concrete can be tailored to sufficiently attenuate all types of ionizing radiation, with particular emphasis on gamma-rays and neutrons [5,6].
To effectively attenuate neutrons, materials with high neutron cross-sections are necessary [23,24].These materials include both light and heavy elements.Heavy elements cause significant energy loss for neutrons through inelastic scattering [25].On the other hand, elastic scattering by hydrogen (or moderator elements) results in considerable energy transfer, causing the neutron to slow down towards the thermal region [24].A thermal neutron is a neutron that is not bound within an atomic nucleus and has an energy of around 25 meV.This energy corresponds to the most likely speed at 270 K in the Boltzmann distribution.Thermal neutrons typically have a significantly larger effective neutron absorption cross-section for a given nuclide compared to fast neutrons (which have an energy > 1 eV).Therefore, in shielding design, our goal is to slow down fast neutrons to the thermal region in order to enhance their absorption [26].Neutron-absorbing substances such as gadolinium can then absorb these thermal neutrons.Gadolinium ( 157 Gd, with a natural abundance of 15.7%) possesses the highest thermal (< 0.1 eV) neutron capture cross-section among all known stable isotopes [27].
While the replacement of natural aggregates with heavyweight aggregates in the concrete mix design seems relatively simple technologically, the inclusion of powdered additives (usually as cement filler) is generally linked with noticeable alterations to the cement hydration process.In recent years, there has been a visible interest in the introduction of nanosized admixtures into cement-based composites, as their performance is superior to their micro-sized counterparts [8][9][10]21].It was already confirmed in other radiation shielding materials such as polymers [28,29], glass [30,31] and clays [32][33][34], that nanosized admixtures exhibit superior attenuating properties than their micro-sized counterparts.Primary factors attributable to this behavior are (1) surface-to-volume ratio and (2) quantum confinement effects [35].Similarly, various comparative studies on cement-based materials have already confirmed the superiority of nanosized PbO [22], WO 3 [22] and Bi 2 O 3 [9,10] over coarser particles.
However, based on the polymer [36,37] and glass [38] technology experiences, in order to properly and effectively improve the radiation shielding of material it is necessary to use a combination of particles countering specific radiation sources [39].It can be achieved either by mixing particular materials separately or by producing advanced functional fillers such as W/Gd 2 O 3 structures [40,41] or cobalt-doped titania nanocomposites [42].To date, the latter approach has not been reported in the literature related to concrete technology.At the same time, the incorporation of bulk mixtures of materials is available only in a few works.For instance, Soni et al. [43] combined WC and B 4 C particles in concrete, while Nikbin et al. [20] proposed the introduction of nano-Bi 2 O 3 and nano-WO 3 particles for targeted improvement of radiation shielding.
Although many available state-of-the-art studies evaluated the attenuating properties [9,10,20] solely, low interest was paid to the hydration, rheological and fresh performances of cement-based composites containing nano-admixtures.Knowledge related to these aspects is critical, as the amount of nano-admixtures used for radiation shielding purposes usually exceeds the optimal value for ensuring proper mechanical and durability performance of cement-based composites.It is widely agreed that, typically, the optimal content of nanoparticles lies within the range of 0.5 to 5% mass of binder [44][45][46].However, for radiation shielding applications, the amount of a single nano-admixture frequently ranges between 3 and 10 wt% [8, 10,[12][13][14]17,47], therefore, higher replacement rates of a binder with nanoparticles might lead to undesirable alterations in the hydration process of cementitious systems.
For instance, both Bi 2 O 3 micro-and nanosized particles were found to dramatically retard the hydration process, which was reflected in extremely low mechanical performance values of cement pastes [8,47].Such performances might prevent the use of Bi 2 O 3 nanoparticles in technologies such as repairing mortars, prefabrication technology or additive manufacturing (3D printing), where early setting times and improved early strengths are crucial.Available literature confirms that both nanosized [8,47] and microsized [48][49][50] Bi 2 O 3 particles substantially retard the hydration process of cement, with a higher magnitude when nanoparticles are present in the cementitious system.A calorimetric study performed by Cendrowski et al. [47] reported that replacement of cement with 5 wt% and 10 wt% of Bi 2 O 3 nanoparticles resulted in delayed occurrences of the second exothermic peak by 73% and 132% when compared to plain cementitious system with a water-to-binder (w/b) ratio equal to 0.4.
Similarly, knowledge of the performance of Gd 2 O 3 -modified cementitious systems is extremely limited to very few research works that did not deal directly with the radiation shielding performance of the material.In contrast to Bi 2 O 3 , previous studies strongly suggested that Gd 2 O 3 particles exhibit a relatively neutral and limited effect on the fresh and hardened properties of cement-based composites [47,51].

Research Significance
To date, there has been a very limited amount of work that has reported attempts to improve both the gamma-ray and neutron shielding performance of cement-based composites by the inclusion of micro-and nanosized additives.In addition, these works were limited only to the incorporation of bulk mixtures of materials.This study aims to fill the gap in the state-of-art and propose, for the very first time, a synthesis method of functional additives for simultaneous improvement of gamma-ray and neutron attenuation in cementbased composites.The proposed structures are composed of Bi 2 O 3 and Gd 2 O 3 particles and SiO 2 coating for addressing both radiation shielding performance and reducing the negative effect of the early setting of cementitious composite.To mitigate the obstacles related to the highly delayed hydration process and strength gain of a cementitious system containing Bi 2 O 3 nanoparticles, a sol-gel coating method was introduced to produce a functional additive.The sol-gel coating is the most popular method for improving the reactivity of nanoparticles.The synthesis of thin reactive silica layers on the surface of particles enhances reactivity and incorporates nanosized admixtures in the cementitious matrix.Such technology was already applied to produce silica-coated Fe 3 O 4 [52], TiO 2 [53], carbon nanotubes [54] as well as steel fibers [55] for cement-based concrete applications.
In this study, two types of silica coatings were proposed.By recommending two synthesis methods, Bi 2 O 3 -Gd 2 O 3 -SiO 2 structures are synthesized with varied silica morphology, reactivity and specific surface area.Afterward, synthesized functional additives are evaluated in standard Portland cement to assess the performance of structures in a cementitious system.Both hydration/early age performance and radiation shielding properties of cement pastes were evaluated.For comparison purposes, control cement paste and cement containing pristine materials (Bi 2 O 3 and Gd 2 O 3 ) were produced.As an outcome, synthesis method of radiation shielding additives for cement mortars and concrete applications was proposed.

Materials
To produce cement pastes, cement CEM I 42.5 R (Góra żd że, Poland) and tap water were used.Commercially available bismuth oxide (Bi 2 O 3 -cat.No. 637017) and gadolinium oxide (Gd 2 O 3 -cat.No. 278513) were purchased from Merck, Poland.The densities of Bi 2 O 3 and Gd 2 O 3 determined using a helium pycnometer (Ultrapyc 1200e, Quantachrome Instruments, Anton Paar Group AB, Boynton Beach, FL, USA) were 8.68 g/cm 3 and 7.28 g/cm 3 , respectively.To produce the silica coatings, tetraethyl orthosilicate (TEOS) was used (Merck, Warsaw, Poland) as a source of silica.Ethylene alcohol (ETOH) and 26% ammonium solution (NH 4 OH) were purchased from Warchem, Poland.Figure 1 presents the particle size distribution determined by the laser diffraction method (Mastersizer2000, Malvern Pananalytical, Malvern, Worcs, UK) of pristine Bi 2 O 3 and Gd 2 O 3 particles and cement.Both powders were found to be finer than cement, with the highest fineness of Bi 2 O 3 .The detailed description of primary particles will be discussed in detail in Section 4.1.
To produce cement pastes, cement CEM I 42.5 R (Górażdże, Poland) and tap w were used.Commercially available bismuth oxide (Bi2O3-cat.No. 637017) and gad ium oxide (Gd2O3-cat.No. 278513) were purchased from Merck, Poland.The dens of Bi2O3 and Gd2O3 determined using a helium pycnometer (Ultrapyc 12 Quantachrome Instruments, Anton Paar Group AB, Boynton Beach, FL, USA) were g/cm 3 and 7.28 g/cm 3 , respectively.To produce the silica coatings, tetraethyl orthosil (TEOS) was used (Merck, Warsaw, Poland) as a source of silica.Ethylene alcohol (ET and 26% ammonium solution (NH4OH) were purchased from Warchem, Poland.Fi 1 presents the particle size distribution determined by the laser diffraction method (M tersizer2000,Malvern Pananalytical, Malvern, Worcs, UK) of pristine Bi2O3 and G particles and cement.Both powders were found to be finer than cement, with the hig fineness of Bi2O3.The detailed description of primary particles will be discussed in d in Section 4.1.

Synthesis Process of Functional Additives
The synthesis of silica coatings was performed using two methods marked as A B. Details on silica synthesis methods were described comprehensively in the auth previous work [47] related to silica coating of singular particles.In both methods, e amounts of the Bi2O3 and Gd2O3 were used.For each singular synthesis silica, regard of method (A or B), a fixed amount of 45 g of Bi2O3 and Gd2O3 (mass ratio 1:1) was us In method A (Figure 2A), 900 mL of ethylene alcohol (ETOH) and 45 mL of am nium solution were mixed using a magnetic stirrer.After 1 h of mixing, the mixtu the Bi2O3 and Gd2O3 was dispersed in the prepared solution.To disperse metal oxid magnetic stirrer and ultrasound were used alternately.Further to the obtained sus sion, 27 mL of TEOS was added into the flask and stirred continuously for 24 h at r temperature.In the ending stage, the suspension was dried in air at 80 °C.
Synthesis of silica type B (Figure 2B) was performed in a mixture of 720 mL of tilled water, 180 mL of ETOH and 45 mL of ammonium solution.A mixture of solv was poured into the flask and stirred using a magnetic stirrer.Further to the stirred ture, 27 mL of TEOS was added.Afterward, the solution obtained was stirred for 2 with the stirrer rotation speed of 500 rpm.After 24 h, metal oxides were added to

Synthesis Process of Functional Additives
The synthesis of silica coatings was performed using two methods marked as A and B. Details on silica synthesis methods were described comprehensively in the authors' previous work [47] related to silica coating of singular particles.In both methods, equal amounts of the Bi 2 O 3 and Gd 2 O 3 were used.For each singular synthesis silica, regardless of method (A or B), a fixed amount of 45 g of Bi 2 O 3 and Gd 2 O 3 (mass ratio 1:1) was used.
In method A (Figure 2A), 900 mL of ethylene alcohol (ETOH) and 45 mL of ammonium solution were mixed using a magnetic stirrer.After 1 h of mixing, the mixture of the Bi 2 O 3 and Gd 2 O 3 was dispersed in the prepared solution.To disperse metal oxides, a magnetic stirrer and ultrasound were used alternately.Further to the obtained suspension, 27 mL of TEOS was added into the flask and stirred continuously for 24 h at room temperature.In the ending stage, the suspension was dried in air at 80 • C.
Synthesis of silica type B (Figure 2B) was performed in a mixture of 720 mL of distilled water, 180 mL of ETOH and 45 mL of ammonium solution.A mixture of solvents was poured into the flask and stirred using a magnetic stirrer.Further to the stirred mixture, 27 mL of TEOS was added.Afterward, the solution obtained was stirred for 24 h, with the stirrer rotation speed of 500 rpm.After 24 h, metal oxides were added to the suspension and stirred for an additional hour.In the last step, the suspension was cast on the flat surface of the metal trays and evaporated in air at 80 • C for 24 h.
suspension and stirred for an additional hour.In the last step, the suspension was cast on the flat surface of the metal trays and evaporated in air at 80 °C for 24 h.

Mixture Design and Specimen Preparation
To assess the performance of synthesized structures, four representative Portland cement paste (Table 1) mixes were designed with a fixed water-to-binder (w/b) ratio of 0.4.In mixtures containing pristine nanoparticles and synthesized structures, a fixed 1:1 mass ratio of Bi2O3 to Gd2O3 was set.In nanomodified mixes, cement was replaced with 10% (by weight) pristine or synthesized structures.Therefore, the content of singular nanoparticles does not exceed a widely agreed state-of-art optimal dosage of 5 wt%.A control cement paste (without additives) specimen was designated as control.
Cement pastes were mixed using an automatic hand mixer.Before mixing, nanoparticles were mixed with tap water (mixing water) and ultrasonicated for 10 min using an ultrasonic bath.The mixing process of cement paste lasted 180 s following the procedure: (i) 30 s slow mixing, (ii) 60 s fast mixing, (iii) 30 s rest, (iv) 60 s fast mixing.After mixing, the material was used to determine its fresh properties, while the rest of the material was used to cast 20 × 20 × 20 mm specimens.After casting, specimens were covered with foil and cured at room temperature for 24 h (T = 20 °C ± 1 °C).Afterward, cubes were demolded and stored in the container up to 28 d with RH = 95% and T = 20 °C ± 1 °C.

Mixture Design and Specimen Preparation
To assess the performance of synthesized structures, four representative Portland cement paste (Table 1) mixes were designed with a fixed water-to-binder (w/b) ratio of 0.4.In mixtures containing pristine nanoparticles and synthesized structures, a fixed 1:1 mass ratio of Bi 2 O 3 to Gd 2 O 3 was set.In nanomodified mixes, cement was replaced with 10% (by weight) pristine or synthesized structures.Therefore, the content of singular nanoparticles does not exceed a widely agreed state-of-art optimal dosage of 5 wt%.A control cement paste (without additives) specimen was designated as control.Cement pastes were mixed using an automatic hand mixer.Before mixing, nanoparticles were mixed with tap water (mixing water) and ultrasonicated for 10 min using an ultrasonic bath.The mixing process of cement paste lasted 180 s following the procedure: (i) 30 s slow mixing, (ii) 60 s fast mixing, (iii) 30 s rest, (iv) 60 s fast mixing.After mixing, the material was used to determine its fresh properties, while the rest of the material was used to cast 20 × 20 × 20 mm specimens.After casting, specimens were covered with foil and cured at room temperature for 24 h (T = 20 • C ± 1 • C).Afterward, cubes were demolded and stored in the container up to 28 d with RH = 95% and T = 20

Nanomaterials Characterization
To assess the morphological properties and composition of pristine and synthesized structures, transmission electron microscopy (TEM, Tecnai F30 with a field emission gun operating at 200 kV, Thermo Fisher Scientific, Waltham, MA, USA) with energy dispersive X-ray spectroscopy-EDS (EDX, FEI Corporation, Hillsboro, OR, USA) was used.The crystallinity and phase composition of nanomaterials were characterized using an Aeris diffractometer (Malvern Pananalytical, Malvern, Worcs, UK) with Cu-Kα radiation (λ = 1.544Å).The BET specific surface area was measured using an N 2 adsorption/desorption isotherm (Micromeritics ASAP 2460, Norcross, GA, USA).Ultrapyc 1200e (Quantachrome) helium pycnometer was used to assess the specific gravity of newly synthesized structures.

Rheological Properties
Consistency of cement pastes (spread diameter) was determined using a mini-cone that was 36 mm (top diameter) × 60 mm (bottom diameter) × 60 mm (height).A similar mini-cone geometry was used in other studies for cement paste purposes [47,56].The mini-cone was placed on the flow table and raised vertically.Afterwards, 15 hits of the flow table were applied and the spread diameter was determined in two perpendicular directions, with the mean value taken as a representative.
The rheological properties of cement pastes were measured using an MCR 72 (Anton Paar) compact rheometer with vane geometry.A ribbed cup with volume (v=120 cm 3 ) was used to prevent a wall slip effect.In the first stage of the test, the mixture was pre-sheared for 30 s at a constant shear rate of 100 s −1 and afterward, the shear stress of mixtures was determined for a total of 300 s at 20 points distributed on a logarithmic scale with shear rates ranging from 100 s −1 to 0.1 s −1 .
Various rheological fitting models such as the linear Bingham model (BM), modified-Bingham model (MBM), and Herschel-Bulkley model (HBM) have been employed to address the non-Newtonian and pseudoplastic behavior of cement-based mixtures [57,58].In this study, the HBM was found to be the most accurate in fitting the data with a minimum correlation coefficient (R 2 ) of 0.99697.The yield shear stress (τ 0 ), flow rheology index (n) and consistency coefficient (K) of each mixture were determined using (Equation ( 1)).In the Herschel-Bulkley model, the consistency coefficient serves as a crucial parameter that gauges a material's ability to resist flow after surpassing the yield stress.This coefficient is indicative of the material's internal structure and degree of internal friction, essentially revealing the material's resistance to deformation.The flow rheology index is another important parameter that characterizes the non-Newtonian behavior of a material.It defines the nature of the material's flow after exceeding the yield stress.This parameter depicts the shear-thickening and shear-thinning behavior of the material.Shear-thickening occurs when the n>1 and shear-thinning occurs when the n < 1 [59].

Isothermal Calorimetry
Hydration of cement paste was monitored using a TAM Air 8-channel calorimeter (TA Instruments, New Castle, DE, USA) at 20 • C for 168 h.For the measurement, the 12.6 g of cement paste was placed in a 20 mL PE vial, tightly sealed with a screw cap and inserted into the calorimeter.The vials with sand were used as a reference.

Thermogravimetric Analysis (TGA)
Thermogravimetric analysis (TGA) of cement pastes after 2 d, 7 d and 28 d of curing was performed using TA Instruments 5500 in nitrogen, with a heat rate 10 • C/min, up to 1000 • C. To perform the measurement, the cubic specimens were crushed and immersed in the solvent (isopropanol) to stop the hydration process for 1 h.Afterward, the material was dried using a vacuum dryer (40 • C) and milled to fine powder.The mass loss caused by the dehydroxylation of calcium hydroxide (CH) in a temperature range of around 400 • C to 500 • C was determined with the tangential method, according to Scrivener et al. [60].
The following equation was used to determine CH content in the paste (Equation ( 2)): where is the mass of calcium hydroxide in the paste, M 400 • C−500 • C is the range of mass loss caused by the dehydroxylation of calcium hydroxide (determined using the tangential method), M 1000 • C mass of specimen at 900 represent the molar masses of calcium hydroxide (74 g/mol) and water (18 g/mol).

Mechanical Performance
The compressive strength of specimens was determined on 20 mm × 20 mm × 20 mm specimens using electromechanical universal tester UNIFRAME 250 (CONTROLS S.p.A., Liscate, Italy) with a loading rate of 0.6 MPa/s.The measurement was performed after 1 d, 2 d, 7 d and 28 d of curing.Three specimens were tested and the mean value was taken as a representative.

Mercury Intrusion Porosimetry
To analyze the microstructural properties, mercury intrusion porosimetry tests were performed using Quantachrome Poremaster 60.The specimens in question (0.7 cm × 0.7 cm × 2.0 cm) were cut from the middle section of cubic samples.For each composite, two tests were conducted to verify the repeatability of the results.The tests were carried out after 28 days of composite curing.Before the tests, the specimens were dried at 70 • C for 48 h.The mercury surface tension was set at 0.48 N/m and the contact angle was set to 140 • upon intrusion.The specimens were placed into measurement cells and filled with mercury in a low-pressure chamber (to 0.34 MPa).Furthermore, the cells were inserted into a pressure chamber and subjected to high pressures (up to 413 MPa).Based on the tests, cumulative and log differential graphs of porosity distribution were made and the basic properties of composites were given, such as porosity, specific surface area and volume density.

Radiation Shielding Performance Test
In this study, two silica-coated Bi 2 O 3 /Gd 2 O 3 structures, A and B, were utilized to produce silica nanospheres.These methods yielded distinct structures while preserving the same theoretical elemental composition.It is theoretically unfeasible to determine the expected difference in radiation capacity due to the fabrication method.Therefore, only the experimental results were considered.The radiation shielding properties of the prepared samples were assessed using narrow-beam transmission geometry.
A comprehensive explanation of the experimental setup was already presented in the previous work of our group [8].In the gamma-ray experiment, we used gammaray lines emitted by radioactive point sources, specifically 133 Ba, 137 Cs, 60 Co and 233 Th.These sources covered an energy range of 81-2614 keV and were used for both the sample irradiation and detector calibration.To record and analyze the obtained spectra, we utilized a 3 × 3 NaI(Tl) scintillator detector, a multichannel analyzer and Genie 2000 software (Canberra).The gamma-ray setup measurements are presented in Figure 3A, which depicts the essential geometry.Furthermore, an investigation on neutrons was performed, utilizing an Am-Be neutron source with an activity of 3.7 GBq.This study was conducted using the experimental setup presented in Figure 3B.
which depicts the essential geometry.Furthermore, an investigation on neutrons was performed, utilizing an Am-Be neutron source with an activity of 3.7 GBq.This study was conducted using the experimental setup presented in Figure 3B.From this geometric configuration and Lambert-Beer law, the following fundamental attenuation coefficients were obtained: the linear coefficient, LAC (µ); mass attenuation coefficient, MAC (µm); total macroscopic cross-section of the slow neutron attenuation coefficient (Σs); half-value layer (HVL); and mean free path (MFP) using the equations listed in Table 2.

Nanomaterials Characterization
The TEM images of gadolinium and bismuth oxide structures after silica coating (Figure 4A,D) show agglomerates made from pristine metal oxide particles.In the agglomerates, two types of particles can be distinguished.Small particles' diameter ranges from 93 nm to 643 nm, with a mean diameter of about 250 nm corresponding to the size of bismuth oxide particles (according to the TEM analysis, Figure 5A,B).According to the TEM images, bigger and darker particles correspond to the size and shape of pristine gadolinium oxide particles.The Gd2O3 crystals' size ranges from 101 nm to 248 nm, with a more spherical and regular shape.The Gd2O3 crystals are in the form of agglomerates From this geometric configuration and Lambert-Beer law, the following fundamental attenuation coefficients were obtained: the linear coefficient, LAC (µ); mass attenuation coefficient, MAC (µ m ); total macroscopic cross-section of the slow neutron attenuation coefficient (Σ s ); half-value layer (HVL); and mean free path (MFP) using the equations listed in Table 2.

Nanomaterials Characterization
The TEM images of gadolinium and bismuth oxide structures after silica coating (Figure 4A,D) show agglomerates made from pristine metal oxide particles.In the agglomerates, two types of particles can be distinguished.Small particles' diameter ranges from 93 nm to 643 nm, with a mean diameter of about 250 nm corresponding to the size of bismuth oxide particles (according to the TEM analysis, Figure 5A,B).According to the TEM images, bigger and darker particles correspond to the size and shape of pristine gadolinium oxide particles.The Gd 2 O 3 crystals' size ranges from 101 nm to 248 nm, with a more spherical and regular shape.The Gd 2 O 3 crystals are in the form of agglomerates with sizes measured between 256 nm and 2000 nm (Figure 5C,D).The XRD analysis of the bismuth and gadolinium oxides after silica enrichment is presented in Figure 4G.No significant difference between the samples enriched with the silica (type A and B) was noticed.All peaks in the samples were assigned to the β-Bi 2 O 3 and Gd 2 O 3 phase, according to the PDF 27-0050 and JCPDS card no.43-1014, respectively.Additionally, the spectra of pristine metal oxides are presented in Figure 4G.The XRD analyses confirm that the synthesis of silica coatings did not influence the chemical structure of the metal oxide core.
Nanomaterials 2024, 14, x FOR PEER REVIEW 9 of 25 with sizes measured between 256 nm and 2000 nm (Figure 5C,D).The XRD analysis o the bismuth and gadolinium oxides after silica enrichment is presented in Figure 4G.No significant difference between the samples enriched with the silica (type A and B) was noticed.All peaks in the samples were assigned to the β-Bi2O3 and Gd2O3 phase, accord ing to the PDF 27-0050 and JCPDS card no.43-1014, respectively.Additionally, the spec tra of pristine metal oxides are presented in Figure 4G.The XRD analyses confirm tha the synthesis of silica coatings did not influence the chemical structure of the metal oxide core.TEM images revealed significant differences between samples enriched with silica according to methods A and B. Silica synthesized employing method A has a solid structure and smooth surface that is free from impurities.The thickness of the silica shell ranges from 26 nm to 35 nm, with a mean thickness of around 30 nm.The TEM image analysis estimated the thickness of the silica shell.High-magnification TEM images of the silica synthesized according to method A are presented in Figure S1A,B (Supplementary Materials) and the previous publication [47].The introduction of silica coatings resulted in a decrement of the structures' specific gravities compared to the values reported for singular particles.It is attributed to the substantially lower density of SiO2 compared to Bi2O3 and Gd2O3 particles (Table 3).
The synthesized silica particles from method B agglomerate at sizes below 15 nm, resembling silica particles (Figure S1 C,D, Supplementary Materials).Such a structure should be characterized by a higher surface area than silica synthesized via method A. Also, in the sample synthesized by method B, additional free silica can be noticed.Other higher magnification TEM images of silica coating over metal oxides are presented in Figure S2 (Supplementary Materials).
The pores and surface area analysis proved the changes in the silica structures.According to BET analysis, Bi2O3 + Gd2O3 with silica synthesized using methods A and B have pore volumes of 0.0023 cm 3 /g and 0.04682 cm 3 /g, respectively.Hence, the surface TEM images revealed significant differences between samples enriched with silica according to methods A and B. Silica synthesized employing method A has a solid structure and smooth surface that is free from impurities.The thickness of the silica shell ranges from 26 nm to 35 nm, with a mean thickness of around 30 nm.The TEM image analysis estimated the thickness of the silica shell.High-magnification TEM images of the silica synthesized according to method A are presented in Figure S1A,B (Supplementary Materials) and the previous publication [47].The introduction of silica coatings resulted in a decrement of the structures' specific gravities compared to the values reported for singular particles.It is attributed to the substantially lower density of SiO 2 compared to Bi 2 O 3 and Gd 2 O 3 particles (Table 3).The synthesized silica particles from method B agglomerate at sizes below 15 nm, resembling silica particles (Figure S1 C,D, Supplementary Materials).Such a structure should be characterized by a higher surface area than silica synthesized via method A. Also, in the sample synthesized by method B, additional free silica can be noticed.Other higher magnification TEM images of silica coating over metal oxides are presented in Figure S2 (Supplementary Materials).
The pores and surface area analysis proved the changes in the silica structures.According to BET analysis, Bi 2 O 3 + Gd 2 O 3 with silica synthesized using methods A and B have pore volumes of 0.0023 cm 3 /g and 0.04682 cm 3 /g, respectively.Hence, the surface area of the samples reaches 8.0129 m 2 /g and 45.0256 m 2 /g, respectively (Table 3).The difference in the total volume of pores and surface area proves the silica coating's microscopic analysis.
Upon comparison of methods A and B, significant differences were observed in the synthesized silica samples under optical microscopic analysis.Silica synthesized via method A has a white color and powder form (Figure 6A,B).This microscopic analysis is in accordance with the previous studies on the synthesis of silica nanomaterials (spheres and coatings) [64][65][66].These studies indicated that silica synthesized via method A has a non-porous and amorphous character [66].
area of the samples reaches 8.0129 m 2 /g and 45.0256 m 2 /g, respectively (Table 3).The difference in the total volume of pores and surface area proves the silica coating's microscopic analysis.Upon comparison of methods A and B, significant differences were observed in the synthesized silica samples under optical microscopic analysis.Silica synthesized via method A has a white color and powder form (Figure 6A,B).This microscopic analysis is in accordance with the previous studies on the synthesis of silica nanomaterials (spheres and coatings) [64][65][66].These studies indicated that silica synthesized via method A has a non-porous and amorphous character [66].On the contrary, silica synthesized using method B is transparent and has a crystal form resembling glass.Figure 6C,D shows optical microscope images of dried silica.The translucent pieces of silica can be noticed as white spots.This white powder is microsized silica from the crashed, broken or grounded silica pieces.Figure 6D shows the image of the white, small particles at higher magnification.This image proves that white particles have a similar structure/form to the macro-sized silica pieces.
The difference between both methods can also be observed within the synthesized silica particles.The solution from the silica synthesized according to method A changes from transparent and colorless to milky white, whereas the solution remains transparent and colorless during the second method.Our analysis was proven by Yun et al. [67], which shows dependence between the size of the particles and the color of the silica On the contrary, silica synthesized using method B is transparent and has a crystal form resembling glass.Figure 6C,D shows optical microscope images of dried silica.The translucent pieces of silica can be noticed as white spots.This white powder is micro-sized silica from the crashed, broken or grounded silica pieces.Figure 6D shows the image of the white, small particles at higher magnification.This image proves that white particles have a similar structure/form to the macro-sized silica pieces.
The difference between both methods can also be observed within the synthesized silica particles.The solution from the silica synthesized according to method A changes from transparent and colorless to milky white, whereas the solution remains transparent and colorless during the second method.Our analysis was proven by Yun et al. [67], which shows dependence between the size of the particles and the color of the silica suspensions.As reported, silica up to 30 nm shows a tendency to create transparent solutions.Further increase to the size of colloidal silica changes color to white and milky white.Their studies are in accordance with our present and previous studies on nanosilica.
The most crucial advantage of method B over method A is the requirements and complexity of synthesis.Synthesis of SiO 2 via method B allows for the reduction of 80% of the ethanol by replacing it with distilled water [47], in contrary to method A [64,65].It allows for a significant decrease in synthesis costs.For synthesizing silica shells, method B offers a considerable advantage over other methods due to its simple and uncomplicated synthesis procedure described in detail in [47].The complexity of method A compared to method B lies in the synthesis step where metal oxides are introduced and additional treatment of metal oxides is required.Synthesized silica shell via method A requires proper homogenization of metal oxides using stirring and sonication in the solution before TEOS hydrolyses and condensation.In method B, metal oxides are introduced after TEOS hydrolyses.

Fresh Properties
The rheology measurements were performed on all the mixtures, and the principal rheology parameters, including yield shear stress, consistency coefficient and rheological index of the mixtures, were calculated using the HBM fitting model (Figure 7).Table 4 shows that all the mixtures have n < 1, suggesting that the mixtures exhibit shear-thinning behavior.
registered for the mixture containing plain BG.The reason has to be associated with the increased interaction between the particles within the mixture.The utilization of method B to coat BG particles with silica powder could introduce additional interparticle forces, such as Van der Waals forces or electrostatic interactions.These increased interactions enhance particle interlocking, resulting in higher yield shear stress.At the same time, the improved particle-particle interaction can increase internal friction and resistance to flow, leading to an increase in consistency coefficient.
Figure 7B depicts a mini-slump flow test measuring fresh concrete's overall consistency or workability.The results indicated that the inclusion of BG particles led to a slight reduction in the spread diameter from 139 mm for the control mix to 133 mm for the BG.
Employing silica coating method A on BG particles resulted in a decrease in spread diameter to 128 mm when compared to the control sample.The highest reduction in spread diameter was registered for the sample coated with method B, with a 14% spread diameter reduction compared to that of the BG mixture.The observation is in line with the previous study [68], confirming that due to large specific surface area of nanosilica the consistency of cementitious composites have a tendency to decrease.The results obtained through mini-slump tests do not correspond with the trends found on dynamic yield shear stress values obtained in the rheology measurements.Previous research [69] has also suggested that basic assessments such as flow table and slump may not accurately read rheological characteristics.

Isothermal Calorimetry
Figure 8 presents a heat flow, and cumulative heat curves of cement pastes tested up to seven days of hydration.In this study, the isothermal calorimetry tests were conducted to investigate the hydration behavior of the mixtures, aiming to identify the reason for the early age compressive strength reduction (Section 4.5).For better readability and discussion, the heat flow results (Figure 8A) were presented up to 72 h.The results indicated that the incorporation of BG additives without coating results in a substantial increase in the yield shear stress of the mixture.As can be seen in Table 4, the yield shear stress of the paste increased from 15.9 Pa for the control sample and reached 53.2 Pa for the BG mixture.The incorporation of Bi 2 O 3 and Gd 2 O 3 nanoparticles in pastes increases yield shear stress.This effect can be attributed to the improved particle packing and interlocking within the paste matrix, which leads to enhanced resistance against deformation and flow.The nanoparticles improve the interparticle forces, resulting in stronger particle-particle interactions and a more compact microstructure.These changes contribute to the increased yield shear stress observed in the pastes.
The coating implementation, however, results in a distinct alternation in the rheological parameters of the fresh mortar.The results indicated that the mixture containing BG particles modified through method A exhibited a reduction in the yield shear stress compared to plain BG additive incorporation, while the consistency coefficient slightly increased.The results indicate that the BG-A showed a decrease in yield shear stress (i.e., about 19%) and a slight increase in consistency coefficient (i.e., 4%) compared to the BG.The coating process could have influenced the dispersion and aggregation behavior of the BG particles.Improved dispersion through method A might have reduced particle aggregation, resulting in a more dispersed state within the paste.Moreover, as confirmed in the TEM results, coating method A produces particles with smooth surfaces free from impurities.Enhanced dispersion coupled with the presence of particles with smooth surfaces decrease the internal friction and facilitate the material flow under shear force, reducing yield shear stress.Moreover, silica coating can act as a lubricant or flow-enhancing additive due to its low friction and smooth surface properties.The presence of silica powder as a coating on BG particles can contribute to reduced friction and resistance to flow, leading to a decrease in yield shear stress.
Coating method B, on the other hand, results in a slight increase in the yield shear stress and consistency coefficient by 6% and 0.5%, respectively, compared to the values registered for the mixture containing plain BG.The reason has to be associated with the increased interaction between the particles within the mixture.The utilization of method B to coat BG particles with silica powder could introduce additional interparticle forces, such as Van der Waals forces or electrostatic interactions.These increased interactions enhance particle interlocking, resulting in higher yield shear stress.At the same time, the improved particle-particle interaction can increase internal friction and resistance to flow, leading to an increase in consistency coefficient.
Figure 7B depicts a mini-slump flow test measuring fresh concrete's overall consistency or workability.The results indicated that the inclusion of BG particles led to a slight reduction in the spread diameter from 139 mm for the control mix to 133 mm for the BG.
Employing silica coating method A on BG particles resulted in a decrease in spread diameter to 128 mm when compared to the control sample.The highest reduction in spread diameter was registered for the sample coated with method B, with a 14% spread diameter reduction compared to that of the BG mixture.The observation is in line with the previous study [68], confirming that due to large specific surface area of nanosilica the consistency of cementitious composites have a tendency to decrease.The results obtained through mini-slump tests do not correspond with the trends found on dynamic yield shear stress values obtained in the rheology measurements.Previous research [69] has also suggested that basic assessments such as flow table and slump may not accurately read rheological characteristics.

Isothermal Calorimetry
Figure 8 presents a heat flow, and cumulative heat curves of cement pastes tested up to seven days of hydration.In this study, the isothermal calorimetry tests were conducted to investigate the hydration behavior of the mixtures, aiming to identify the reason for the early age compressive strength reduction (Section 4.5).For better readability and discussion, the heat flow results (Figure 8A) were presented up to 72 h.
A remarkable delay in the hydration process of BG paste, when compared to control paste, can be noticed.The occurrence of the second hydration peak in BG was delayed by 71%, and the magnitude of the peak was decreased by 35% when compared to the Control, as shown in Table 5.Similarly, the cumulative heat release of the BG specimen was found to be substantially decreased, especially in the first 48 h of hydration.When compared to the control specimen, the BG paste released only 67%, 55% and 61% percentage of heat emitted by the control paste after 12 h, 16 h and 24 h, respectively.After 168 h, the total heat released from the BG specimen was 8% lower compared to the control paste, which can be attributed to the mass replacement rate of cement with BG.The primary factor responsible for the delayed hydration of cementitious systems containing BG is related to the retarding effect of Bi 2 O 3 nanoparticles.Cendrowski et al. [47] confirmed that a cementitious system comprising 5 wt% and 10 wt% cement replacement with nanosized Bi 2 O 3 delayed the occurrence and power of a second exothermic peak from 12 h 54 min (power of 2.37 mW/g −1 ) to 22 h 21 min (power of 1.65 mW/g −1 ) and 29 h 53 min (power of 1.48 mW/g −1 ), respectively.A remarkable delay in the hydration process of BG paste, when compared to control paste, can be noticed.The occurrence of the second hydration peak in BG was delayed by 71%, and the magnitude of the peak was decreased by 35% when compared to the Control, as shown in Table 5.Similarly, the cumulative heat release of the BG specimen was found to be substantially decreased, especially in the first 48 h of hydration.When compared to the control specimen, the BG paste released only 67%, 55% and 61% percentage of heat emitted by the control paste after 12 h, 16 h and 24 h, respectively.After 168 h, the total heat released from the BG specimen was 8% lower compared to the control paste, which can be attributed to the mass replacement rate of cement with BG.The primary factor responsible for the delayed hydration of cementitious systems containing BG is related to the retarding effect of Bi2O3 nanoparticles.Cendrowski et al. [47] confirmed that a cementitious system comprising 5 wt% and 10 wt% cement replacement with nanosized Bi2O3 delayed the occurrence and power of a second exothermic peak from 12 h 54 min (power of 2.37 mW/g −1 ) to 22 h 21 min (power of 1.65 mW/g −1 ) and 29 h 53 min (power of 1.48 mW/g −1 ), respectively.Contrary to Bi2O3, available literature suggests that lower dosages of Gd2O3 have either a neutral or slightly accelerating effect on the cement hydration process.In comparison, higher dosages (up to 10 wt%) could lead to slight retardation of the cement hydration process [47].Therefore, based on our study, it can be concluded that the primary cause of retardation of the hydration process of BG cement paste is the presence of Bi2O3 in the cementitious system.
The incorporation of silica-coated Bi2O3-Gd2O3 structures resulted in a remarkable acceleration of the hydration process of cementitious systems when compared to BG and, thus, a substantial reduction of the negative impact of Bi2O3 presence in the paste (Figure 9).Compared to the control paste, the occurrence of a second exothermic peak was delayed only by 15% and 14%.In comparison, the peak's magnitude decreased by 17% and 21% for specimens BG-A and BG-B, respectively (Table 5).Contrary to the BG  Contrary to Bi 2 O 3 , available literature suggests that lower dosages of Gd 2 O 3 have either a neutral or slightly accelerating effect on the cement hydration process.In comparison, higher dosages (up to 10 wt%) could lead to slight retardation of the cement hydration process [47].Therefore, based on our study, it can be concluded that the primary cause of retardation of the hydration process of BG cement paste is the presence of Bi 2 O 3 in the cementitious system.
The incorporation of silica-coated Bi 2 O 3 -Gd 2 O 3 structures resulted in a remarkable acceleration of the hydration process of cementitious systems when compared to BG and, thus, a substantial reduction of the negative impact of Bi 2 O 3 presence in the paste (Figure 9).Compared to the control paste, the occurrence of a second exothermic peak was delayed only by 15% and 14%.In comparison, the peak's magnitude decreased by 17% and 21% for specimens BG-A and BG-B, respectively (Table 5).Contrary to the BG specimen, pastes containing synthesized silica-coated structures (BG-A and BG-B), over the entire time of testing, exhibited comparable heat emission to that of the control paste.The emitted heat after 168 h was found to be 3% and 6% lower for specimens BG-A and BG-B, respectively.
A positive effect of newly synthesized structures on the hydration process of cementitious systems can be found, and can be attributed to thin silica-layer coating.Due to the extremely high reactivity of silica, a substantial acceleration of the hydration process is obtained along with improved embedment of the material in the cementitious matrix [70][71][72].Due to the synthesis process and higher amorphicity level of the solid silica shell (type A coating), higher improvement rates were recorded despite a substantially lower surface area of the material compared to a structure containing type B silica coating, which was further confirmed by TG analysis.
the entire time of testing, exhibited comparable heat emission to that of the control paste.The emitted heat after 168 h was found to be 3% and 6% lower for specimens BG-A and BG-B, respectively.A positive effect of newly synthesized structures on the hydration process of cementitious systems can be found, and can be attributed to thin silica-layer coating.Due to the extremely high reactivity of silica, a substantial acceleration of the hydration process is obtained along with improved embedment of the material in the cementitious matrix [70][71][72].Due to the synthesis process and higher amorphicity level of the solid silica shell (type A coating), higher improvement rates were recorded despite a substantially lower surface area of the material compared to a structure containing type B silica coating, which was further confirmed by TG analysis.

Thermogravimetric Analysis (CH Content)
The results of CH content determined from TG curves are presented in Figure 10.Lower CH content in specimens containing additives can be found after 2 d and 7 d due to a more deficient amount of cement available for hydration when compared to the control paste.The inclusion of synthesized structures resulted in lower CH contents in the BG-A and BG-B pastes as a result of the pozzolanic activity of the silica shell.It is widely agreed that nanosized silica exhibits spectacular pozzolanic activity, and this effect is highly related to various parameters including amorphicity, the size of the particle and its specific surface area [73][74][75][76].The lowest CH content was reported for BG-A, followed by BG-B, which is in line with the activity of additives in the calorimetric study.After 2 d, the CH content in BG and BG-B specimens was comparable, while BG-A exhibited 24% and 7% lower CH content than control and BG, respectively.After 7 d, BG-A and BG-B showed 18% and 13% lower CH content when compared to BG.After 28 days of curing, the discrepancy between CH content in specimens has been increased, confirming the superior pozzolanic activity of BG-A particles attributed to the highly amorphous level of type A silica coating compared to type B coating.

Thermogravimetric Analysis (CH Content)
The results of CH content determined from TG curves are presented in Figure 10.Lower CH content in specimens containing additives can be found after 2 d and 7 d due to a more deficient amount of cement available for hydration when compared to the control paste.The inclusion of synthesized structures resulted in lower CH contents in the BG-A and BG-B pastes as a result of the pozzolanic activity of the silica shell.It is widely agreed that nanosized silica exhibits spectacular pozzolanic activity, and this effect is highly related to various parameters including amorphicity, the size of the particle and its specific surface area [73][74][75][76].The lowest CH content was reported for BG-A, followed by BG-B, which is in line with the activity of additives in the calorimetric study.After 2 d, the CH content in BG and BG-B specimens was comparable, while BG-A exhibited 24% and 7% lower CH content than control and BG, respectively.After 7 d, BG-A and BG-B showed 18% and 13% lower CH content when compared to BG.After 28 days of curing, the discrepancy between CH content in specimens has been increased, confirming the superior pozzolanic activity of BG-A particles attributed to the highly amorphous level of type A silica coating compared to type B coating.

Mechanical Performance
Table 6 presents the compressive strength values after 1 d, 2 d, 7 d and 28 days of curing.For more efficient discussion on early age characteristics of the material, compressive strength values were additionally presented as a percentage of the control (C) specimen's strength at the selected testing date (Figure 11).

Mechanical Performance
Table 6 presents the compressive strength values after 1 d, 2 d, 7 d and 28 days of curing.For more efficient discussion on early age characteristics of the material, compressive strength values were additionally presented as a percentage of the control (C) specimen's strength at the selected testing date (Figure 11).
As seen in Figure 11, for all the compositions containing BG particles, with and without silica coating, the compressive strength is lower than that of the corresponding control sample at early hydration ages (i.e., 1, 2 and 7 days of curing).On the first day of the hydration process, paste mixes BG, BG-A and BG-B attained 17.9%, 53.8% and 45.1% compressive strength of control samples, respectively.As the hydration process progressed, the compressive stress of the mixes mentioned above

Mechanical Performance
Table 6 presents the compressive strength values after 1 d, 2 d, 7 d and 28 days of curing.For more efficient discussion on early age characteristics of the material, compressive strength values were additionally presented as a percentage of the control (C) specimen's strength at the selected testing date (Figure 11).
As seen in Figure 11, for all the compositions containing BG particles, with and without silica coating, the compressive strength is lower than that of the corresponding control sample at early hydration ages (i.e., 1, 2 and 7 days of curing).On the first day of the hydration process, paste mixes BG, BG-A and BG-B attained 17.9%, 53.8% and 45.1% compressive strength of control samples, respectively.As the hydration process progressed, the compressive stress of the mixes mentioned above As seen in Figure 11, for all the compositions containing BG particles, with and without silica coating, the compressive strength is lower than that of the corresponding control sample at early hydration ages (i.e., 1, 2 and 7 days of curing).
On the first day of the hydration process, paste mixes BG, BG-A and BG-B attained 17.9%, 53.8% and 45.1% compressive strength of control samples, respectively.As the hydration process progressed, the compressive stress of the mixes mentioned above gradually improved and approached that of the control sample.After seven days of curing, the compressive strength of the BG, BG-A and BG-B reached 80.5%, 90% and 93.9% compressive strength of the control samples, respectively.Test results revealed that the silica coating approaches, both methods A and B, exhibited a positive effect, diminishing the delay in the hydration process and resulting in a higher compressive strength at early ages.However, method A showed slightly higher mechanical performance, i.e., 19.2%, 5.7% and 5.9% at 1 d, 2 d and 7 days of curing, respectively, than method B. Several authors [77][78][79] have reported that the impact of various nanoparticle incorporation on the properties of cementitious composites is most pronounced during the early stages of hydration.Nevertheless, acquiring adequate mechanical performance is of utmost importance when considering their practical applications [80].According to the study's findings, the inclusion of BG nanoparticles, with or without silica coatings, exhibited no reduction in compressive strength, and all specimens showed comparable value of compressive strength after 28 days of curing despite a decline of cement content in the mix by 10 wt%.

Mercury Intrusion Porosimetry
Paste porosity tests performed using MIP indicate comparable pore distributions in all tested specimens.Comparing the results of porosity distribution (Figure 12) of the tested composites, a slight shift of the pore peak between the control sample and the other pastes can be distinguished (0.075 µm → 0.058 µm).In addition, the BG-B specimen is characterized by the highest value of the pore peak in the range from 0.05 µm to 0.07 µm.In the range from 0.003 µm to 0.02 µm, the pore content in the samples with the addition of metal oxides is clearly higher compared to the control paste, which significantly affected the total pore content and the specific surface area.
findings, the inclusion of BG nanoparticles, with or without silica coatings, exhibited no reduction in compressive strength, and all specimens showed comparable value of compressive strength after 28 days of curing despite a decline of cement content in the mix by 10 wt%.

Mercury Intrusion Porosimetry
Paste porosity tests performed using MIP indicate comparable pore distributions in all tested specimens.Comparing the results of porosity distribution (Figure 12) of the tested composites, a slight shift of the pore peak between the control sample and the other pastes can be distinguished (0.075 µm → 0.058 µm).In addition, the BG-B specimen is characterized by the highest value of the pore peak in the range from 0.05 µm to 0.07 µm.In the range from 0.003 µm to 0.02 µm, the pore content in the samples with the addition of metal oxides is clearly higher compared to the control paste, which significantly affected the total pore content and the specific surface area.Table 7 summarizes the fundamental property values obtained from the MIP measurements.The addition of metal oxides increased the porosity from 1.1% (BG-B) to 2.9% (BG) compared to the control material.It also clearly influenced the value of the specific surface area, which in all cases was higher than in the control sample.The highest surface area was observed in the composite in which the additives were added directly without the silicon coating (BG).The use of a silicon coating in both methods A and B also resulted in an increase in the specific surface area of 21.4% (method A) and 21.8%, respectively (method B).A slightly lower volume density also characterized composites with the addition of oxides compared to the control paste.The obtained results are in line with the theoretical expectations that the replacement of 10 wt% of cement with non-reactive material of higher density will result in higher free water available in the mixture, thus increasing the porosity of the matrix.On the contrary, the inclusion of particles finer than cement resulted in a slight alteration of the pore structure.It shifted towards finer pores in the matrix, resulting in a higher total surface area.In addition, specimens coated with silica resulted in the improvement of the microstructure when compared to BG as a result of higher pozzolanic reactivity of nanosilica, and thus the formation of more refined and compact microstructure, as reported in previous studies [81][82][83].Figure 13A illustrates the variation in the LAC with photon energy for each composition.As the photon energy increased, gamma-ray attenuation decreased.The LAC is influenced by both the incident photon energy and the metal coating method employed on the specimen.On the other hand, the BG sample, which contains pristine nanoparticles of Bi 2 O 3 and Gd 2 O 3 , exhibits improved gamma-ray attenuation compared to the control sample.This improvement is due to the highest contributions of Bi and Gd elements.
ticles finer than cement resulted in a slight alteration of the pore structure.It shifted towards finer pores in the matrix, resulting in a higher total surface area.In addition, specimens coated with silica resulted in the improvement of the microstructure when compared to BG as a result of higher pozzolanic reactivity of nanosilica, and thus the formation of more refined and compact microstructure, as reported in previous studies [81][82][83].Figure 13A illustrates the variation in the LAC with photon energy for each composition.As the photon energy increased, gamma-ray attenuation decreased.The LAC is influenced by both the incident photon energy and the metal coating method employed on the specimen.On the other hand, the BG sample, which contains pristine nanoparticles of Bi2O3 and Gd2O3, exhibits improved gamma-ray attenuation compared to the control sample.This improvement is due to the highest contributions of Bi and Gd elements.Figure 13A also showed that the BG-B sample, containing Bi 2 O 3 and Gd 2 O 3 coated with silica synthesized using method B, demonstrated higher values compared to the BG-A sample.It can be attributed to nanomaterials' larger specific surface area, additional free silica, lower total porosity percentage and higher density compared to silica synthesized using method A, as discussed in Section 4.1 and listed in Table 7.As all of these parameters increase the probability of radiation interaction, the BG-B sample performs more efficiently in gamma-ray attenuating than BG-A.
The improvement values listed in Table 8 were calculated for each sample to quantify the increase in attenuation capacity compared to the control sample.According to Table 8, the addition of pristine nanoparticles of Bi 2 O 3 and Gd 2 O 3 into BG resulted in a significant increase in shielding performance of about 112% at 80 keV.Additionally, both BG-A and BG-B showed a consistent increase in the linear attenuation coefficient compared to the control and BG samples.Furthermore, BG-B demonstrated an approximately 13.5% improvement in shielding performance compared to BG-A, specifically at 80 keV photon energy.The decrease in improvement values as photon energy increases indicates the increasing dominance of Compton scattering interaction with energy.
The half-value layer (HVL) measures the thickness of the specimen required to reduce the intensity of the incident radiation by half.Therefore, it is useful in assessment of the shielding capacity of any material, particularly for engineering design.Figure 13B displays the HVL data, indicating that the BG-B mix had a lower HVL, which demonstrates its enhanced shielding effectiveness.The gap between the curves in Figure 13 shows a measurable improvement in the shielding capacity.

Slow Neutron Attenuation Performance
Table 9 presents the slow neutron linear attenuation coefficients.The inclusion of pristine Bi 2 O 3 /Gd 2 O 3 nanoparticles in the cement paste (BG sample) led to notable improvements in the neutron attenuation capacity, which is consistent with similar findings in gamma-ray studies.Table 9 indicates that the cement pastes with Bi 2 O 3 /Gd 2 O 3 additives exhibited a more significant linear neutron attenuation coefficient (Σ s ) than the control sample.This increase can be attributed to the contribution of Gadolinium (Gd), which possesses a high thermal neutron capture cross-section and increases the density.Furthermore, Table 9 reveals that the choice of silica-coated Bi 2 O 3 /Gd 2 O 3 structures resulted in varying levels of improvement in shielding performance.Similar to the gamma-ray studies, silica-coated structure B exhibited a more significant increase in shielding capacity than silica-coated structure A. Using an approach similar to what is shown in Table 8, the improvement values in Table 9 were computed for each mix and compared with the adjacent sample to determine the increase in attenuation capability.Specimens BG, BG-A and BG-B consistently showed higher Σ s values compared to samples BG, BG-A and BG-B.Specifically, method B demonstrated an approximately 34.8% enhancement in slow neutron shielding when compared to coating method A, as indicated in Table 9.Similar enhancements were also observed in terms of HVL and MFP, suggesting an overall increase in slow neutron shielding.Thus, introducing silica-coatings can ultimately enhance the neutron and gamma shielding capabilities, resulting in improved radiation shielding performance.

Conclusions
Based on the presented research, the following conclusions can be drawn from this work: 1.
Two types of silica-coated Bi 2 O 3 /Gd 2 O 3 structures were synthesized and varied in terms of coating structures, thickness, porosity and surface area, allowing the modification of selected properties such as hydration process, rheology, early strength development and radiation shielding.Proposed structures were found to be beneficial both to early age strength development and radiation shielding properties of cement pastes.Thanks to the introduction of high-Z materials Bi 2 O 3 /Gd 2 O 3 along with the material of neutron capture cross-section (Gd 2 O 3 ), simultaneous improvement of both gamma-ray and neutron attenuation can be achieved.2.
Higher reactivity of synthesized structures, and thus higher acceleration of hydration process and strength development, was obtained in the case of type A silica coating of Bi 2 O 3 /Gd 2 O 3 .Due to the pozzolanic activity of silica coating, lower CH contents were found in specimens containing silica-coated structures when compared to specimens containing uncoated particles.Moreover, higher microstructure refinement in specimens BG-A and BG-B was found compared to specimen BG.Although method B shows slightly lower reactivity in the cement hydration process, it offers advantages related to the additional alteration of rheological properties (increased yield shear stress) which can be useful in selected applications such as 3D printing.At the same time, coating type B has the advantage of lower synthesis time and cost reduction (reduction of ethanol content required for synthesis by 80%).

3.
Rheology tests are an essential and effective tool for selecting nanoparticle candidates and determining the most suitable silica coating methods for specific applications.The incorporation of Bi 2 O 3 and Gd 2 O 3 particles in pastes increases yield shear stress and consistency coefficient.By employing the silica coating method A on BG particles, the yield shear stress of the resulting pastes decreased by around 19%.On the contrary, the application of method B resulted in a slight increase of 6% compared to the plain BG mixture.These findings highlight the significant impact of silica coating on the properties of BG particles, which can ultimately enhance the performance of resulting composites in the hardened state.4.
Both pristine and silica-coated structures can be effectively used as a cement filler at a rate of 10 wt% without deterioration of 28 days compressive strength.However, due to the extremely retarding effect of the hydration process attributed to pristine Bi 2 O 3 , early age performance of cement pastes containing uncoated structures is deficient.Therefore, the introduction of silica coatings can overcome the retardation of hydration, leading to higher hydration heat and strength gain, especially in the first two days of hydration.Thus, specimens BG-A and BG-B exhibited 300% and 251% (after one day) and 25% and 18% (after two days) higher compressive strength than BG specimens.5.
The incorporation of pristine Bi 2 O 3 and Gd 2 O 3 to cement pastes enhanced radiation shielding as expected.The introduction of silica-coated structures resulted in further improvement of the shielding performance of specimens.However, the silica-coated structure synthesized by method B is superior for radiation shielding.The specimens BG-A and BG-B exhibited 7.5% and 13.5% higher LAC, respectively, than the BG specimen at a photon energy of 80 keV.Additionally, the performance of slow neutrons was successfully improved by silica-coated structures.For instance, specimens BG-A and BG-B showed 7.9% and 34.8% higher Σ s , respectively, compared to the BG specimen.6.
As an outcome, silica-coated structures can be successfully used in cement-based composites with demanding early age performances, e.g.repairing mortars, prefabrication technology or additive manufacturing (3D printing) where early setting times and high early strengths are crucial factors for these technologies.

Figure 2 .
Figure 2. Scheme of the particle coating with the silica according to method A (A) and method B (B).

Figure 2 .
Figure 2. Scheme of the particle coating with the silica according to method A (A) and method B (B).

Figure 3 .
Figure 3. Experimental setup for gamma-ray (A) and slow neutron (B) transmission.Reproduced from Sikora et al. [8] with permission from Elsevier, 2021

Figure 4 .
Figure 4. TEM and STEM images of gadolinium and bismuth oxide covered with silica synthesized via methods A (A-C) and B (D-F).The white bar on the STEM images correspond to the length of 500 nm.XRD spectrum (G) shows signals from the gadolinium and bismuth oxide coated with the silica (via different methods) and pristine metal oxides.

Figure 4 .
Figure 4. TEM and STEM images of gadolinium and bismuth oxide covered with silica synthesized via methods A (A-C) and B (D-F).The white bar on the STEM images correspond to the length of 500 nm.XRD spectrum (G) shows signals from the gadolinium and bismuth oxide coated with the silica (via different methods) and pristine metal oxides.

Figure 6 .
Figure 6.Microscopic images of the silica synthesized via method A (A,B) and method B (C,D).The red scale bar in images (A-D) represents 1 mm.

Figure 6 .
Figure 6.Microscopic images of the silica synthesized via method A (A,B) and method B (C,D).The red scale bar in images (A-D) represents 1 mm.

Figure 7 .
Figure 7. Rheology test results (A) and consistency (mean diameter) determined with a mini cone using the flow table method (B).

Figure 7 .
Figure 7. Rheology test results (A) and consistency (mean diameter) determined with a mini cone using the flow table method (B).

Nanomaterials 2024 , 25 Figure 8 .
Figure 8.Heat flow (A) and cumulative heat (B) released from pastes up to 168 h of hydration.

Figure 8 .
Figure 8.Heat flow (A) and cumulative heat (B) released from pastes up to 168 h of hydration.

Figure 9 .
Figure 9. Percentage of cumulative heat released from pastes at different hydration times compared to control paste (C).

Figure 9 .
Figure 9. Percentage of cumulative heat released from pastes at different hydration times compared to control paste (C).

Figure 11 .
Figure 11.Compressive strength of modified specimens in percentage in relation to the control sample.

Figure 11 .
Figure 11.Compressive strength of modified specimens in percentage in relation to the control sample.

Figure 11 .
Figure 11.Compressive strength of modified specimens in percentage in relation to the control sample.

Figure 12 .
Figure 12.Cumulative (A) and log-differential (B) pore size distributions of the tested composites.Figure 12. Cumulative (A) and log-differential (B) pore size distributions of the tested composites.

Figure 12 .
Figure 12.Cumulative (A) and log-differential (B) pore size distributions of the tested composites.Figure 12. Cumulative (A) and log-differential (B) pore size distributions of the tested composites.

Figure 13 .
Figure 13.Gamma-ray attenuation parameters of cement paste as a function of photon energy: LAC (A) and HVL (B).

Table 1 .
The ratio of the cement paste components (wt%).

Table 3 .
Specific gravity and properties of synthesized structures determined from BET.

Table 3 .
Specific gravity and properties of synthesized structures determined from BET.

Table 5 .
Maximum exothermic heat rate and peak occurrence.

Table 5 .
Maximum exothermic heat rate and peak occurrence.

Table 6 .
Compressive strength evolution of cement pastes.

Table 6 .
Compressive strength evolution of cement pastes

Table 7 .
The basic properties of the composited acquired from MIP.

Table 7 .
The basic properties of the composited acquired from MIP.

Table 8 .
Experimental LAC for specimens' control, BG, BG-A and BG-B.LAC improvement factor IMF% of specimens BG, BG-A and BG-B relative to those of samples control, BG, BG-A and BG-B, respectively.Experimental uncertainty ≤ ±3%.